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Abstract 

This report summarizes the initial modeling of the global response of the Bigelow 
Expandable Activity Module (BEAM) to micrometeorite and orbital debris (MMOD) 
impacts using a structural, nonlinear, transient dynamic, finite element code. These 
models complement the on-orbit deployment of the Distributed Impact Detection 
System (DIDS) to support structural health monitoring studies. Two global models 
were developed. The first focused exclusively on impacts on the soft-goods [ fabric - 
envelop) portion of BEAM. The second incorporates the bulkhead to support 
understanding of bulkhead impacts. These models were exercised for random impact 
locations and responses monitored at the on-orbit sensor locations. The report 
concludes with areas for future study. 


1. Introduction 

Micrometeorite and Orbital Debris [MMOD) impacts are a common threat for 
human and robotic spacecraft traveling in low earth orbit, see Ref [1]. Design 
approaches for protecting spacecraft against MMOD impacts that degrade 
performance or cause catastrophic destruction have been studied extensively, see 
Ref. [2]. Flexible spacecraft structural components present a particularly challenging 
structure to protect. A human-rated space habitat demonstration module has been 
fabricated utilizing the design approach documented in a US Patent, see Refs. [3 and 
4]. This demonstration module has been designated the Bigelow Expandable 
Activity Module [BEAM], which will be attached to the International Space Station 
[ISS] as illustrated in Figure 1, see Ref. [5]. 

A number of numerical methods and associated hardware testing for damage 
detection for MMOD type impacts have been documented in the literature, see Ref. 
[6-12]. Historically, hypervelocity testing has been conducted to determine MMOD 
impact performance. Hypervelocity testing is very expensive and is only feasible for 
shields in the latter phases of design. Therefore, advancement of analytical methods 
is important to develop protection against MMOD impacts for novel structures in 
the early trade study and design phases. 

Complementary to advances in soft-goods design and fabrication is progress in the 
development of numerical simulation tools for complex structural systems. For 
example, simulations can incorporate structural aspects such as geometrically 
accurate models and advanced material models that include nonlinear stress-strain 
behaviors, woven fabrics, and inflation. This was demonstrated for the Orion 
Landing System - Advanced Development Project, where the spacecraft landing 
system effectively incorporated modeling of fabric airbags (soft-goods) in the early 
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design process, see Refs. [13 and 14]. Additional simulations incorporating 
structural members with fabrics can be found in Refs. [15-17]. 

The focus here is the global dispersion of responses, due to localized excitations at 
varying locations for structures, where the primary, load-bearing structure is 
composed of tensioned fabric materials. The simulations presented here 
complement the on-orbit deployment of the Distributed Impact Detection System 
(DIDS) structural health monitoring system, see Refs. [18-20], on BEAM. This report 
documents the model development and sample responses for MMOD impact on the 
BEAM structure. First, the impact models will be described. This will be followed by 
numerical results. The report concludes with areas identified for future studies. 

2.0 Numerical models 

Three NASTRAN finite element models were provided by the project, see Ref. [21]. 
These models represented the BEAM structure in various states and were 
designated as "Packed”, "Deflated”, and "Inflated”. The Inflated Model, see Figure 2, 
was selected as the basis for developing the transient models utilized for the results 
reported here. Specifically, two transient-dynamic, finite element models were 
developed to assess the global transient dynamic responses. These models, denoted 
as the "Soft-Goods” and the "Bulkhead" models, will be described in following two 
sections. Both models were executed in LS-DYNA, a commercial, nonlinear, transient 
dynamic finite element simulation tool, see Ref. [22]. 

2.1 Soft-Goods Model 

For the Soft-Goods Model, only the portion of the Inflated Model representing the 
restraint-layer and shield were retained, see Figure 3(a]. The Inflated Model parts 
not incorporated in the Soft-Goods Model are shown in Figure 3[b). The Soft-Goods 
Model contained 12,101 nodes and 11,040 4-node shell elements. For this 
configuration, the nodes attaching the restraint-layer to the Adapter and Passive 
Common Berthing Mechanism were fully restrained in translation. The nodes 
attaching the Soft-Goods to the bulkhead were constrained to behave like a rigid 
body. Concentrated masses were located at the bulkhead interface nodes to equally 
distribute the bulkhead mass. 

The Soft-Goods were assigned fabric material properties. This fabric material model 
implemented in LS-DYNA was originally developed for automotive airbag 
applications. Specifically, the fabric material model is based on an existing 
composite material model, however it is only valid for membrane elements. For 
these simulations, the shell elements were 0.254 cm thick, with nominal material 
properties of elastic modulus=8.96xl0 10 dyne/cm 2 ; shear modulus=3.48x!0 10 
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dyne/cm 2 ; Poisson’s ratio=0.3; and density=11.58 g/cm 3 . These properties were 
derived based on information provided in Ref. [21] as well as private 
communications between the authors and project engineers. To improve numerical 
stability, this material model supports inclusion of a liner that allows for small 
compressive loads. The liner for this application is set to 2% of the restraint layer 
thickness, with the same material properties and a liner damping of 5% critical. 

The total mass of the Soft-Goods Model is 1,468 kg, distributed as 756 kg for the 
restraint layer, 358 kg for the bulkhead, and 354 kg for the adapter and berthing 
mechanism. The inflation pressure is implemented by slowly ramping up the 
pressure on the interior surfaces of the elements over 1 second to a nominal value of 
1.048xl0 6 dyne/cm 2 . Simultaneously, a nodal load is applied in the axial direction to 
the nodes attached to the bulkhead interface to represent the interior pressure 
acting on the bulkhead. Following the inflation load ramp, the impact was 
approximated by a 0.002s triangular force pulse applied to a single node. The 
transient dynamic simulations executed using the Soft-Goods Model were 
completed in less than 10 minutes. 

2.2 Bulkhead Model 

For impacts on the bulkhead, additional parts were incorporated into the Soft-Goods 
Model. The full model and a cross-section are shown in Figure 4. The remainder of 
the bulkhead structure was a direct translation of the NASTRAN model, including 
the linear-elastic material properties assigned. Transmission of the impact energy 
across the Soft-Goods/bulkhead interface was severely attenuated. For that reason, 
the bulkhead model was only used to simulate bulkhead impacts, with attention 
paid to the bulkhead acceleration responses at the DIDS locations also shown on the 
figure. The transient dynamic simulations executed using the Bulkhead Model 
required over 4 hours. 

3.0 Results 

3.1 Modal Comparison of Inflated and Soft-Goods Models 

The global vibration modes of the Soft-Goods Model were evaluated to understand 
the implication of significantly simplifying the Inflated Model mesh. Specifically, a 
comparison of the first three modes for three different models is provided in Figure 
5. The first column represents the modes of the full Inflated Model when executed in 
NASTRAN. The second column represents the modes of the NASTRAN model for the 
Soft-Goods portion only. The third column represents the modes of the Soft-Goods 
LS-DYNA Model. There is less than 3% difference between the NASTRAN and LS- 
DYNA Soft-Goods Models. The Inflated Model modes are 5% lower for the first two 
bending modes with a 26% difference for the third or bouncing mode. This level of 
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agreement was considered sufficient to continue use of the Soft-Goods Model for 
studies of impacts that focus on the restraint layer responses. 

3.2 Soft-Goods Model Studies 

A number of studies were completed using the Soft-Goods Model to further 
understand the implication of some of the modeling unknowns on the acceleration 
responses. The results are presented as resultant acceleration contours with the 
excitation node located at the center. The first parameter studied is the effect of the 
internal pressure, see Figure 6, which directly relates to the restraint layer tension. 
In this case, the baseline is the design condition, with half and double inflation 
pressures also simulated. In the figure, the resultant accelerations are plotted for 3 
times, namely, 0.01, 0.02, and 0.03 seconds after the simulated impact. The wave 
propagation increases with pressure and therefore tension as would be expected. In 
addition, the uniform internal pressure produces faster wave speeds in the 
circumferential direction when compared to the axial direction due to higher fabric 
tension in the circumferential direction. Next, the impact of the elastic modulus on 
the response was studied, see Figure 7. Little effect of the variation in elastic 
modulus was observed on the resultant acceleration. 

3.3 Simulated DIPS responses using Soft-Goods Model 

This section illustrates the method that will be used to generate a response library 
to support on-orbit assessments. A companion report describes the component 
model and contains additional information about through-the-thickness modes. 
Nodes at the locations of the on-orbit DIDS accelerometers on the Soft-Goods Model 
were identified, see Figure 8. Ten random nodes were excited to approximate 
MMOD impacts. Four of the impulse locations were selected for more detailed study, 
see Figure 9. Sample acceleration DIDS time history results for an excitation at each 
location are provided in Figure 10. The Time-of-Arrival of an array of signals can be 
utilized for identification of impact location. The Time-of-Arrival for this 
demonstration was computed as when the acceleration amplitude was greater than 
0.01 g. For the selected excitation locations, the Time-of-Arrival was plotted against 
the geodesic distance from the impulse point for each of the DIDS sensor locations, 
see Figure 11. The geodesic distance was approximated using a cylinder of 127.7 cm 
radius. The cylinder was selected because the geodesic length is a readily- 
programmed closed form expression. When the arrival threshold of the acceleration 
amplitude was increased from 0.01 g, the correlation coefficient of the time of 
arrival versus geodesic distance was substantially less than 1. The average wave 
speed computed from these results is 901 m/s. This process can be expanded to 
generate a large response library. In addition, acceleration history results can be 
used to support development of DIDS impact location identification. 
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3.4 Simulated DIPS responses using Bulkhead Model 

Less is known about soft-goods modeling and propagation of impacts as compared 
to MMOD impacts on solid, metallic structures, see Ref. [2]. Therefore, the focus thus 
far has centered on the simulated soft-goods impacts and responses. Nonetheless, 
sample responses for an impact on the bulkhead have been included. For the case 
illustrated in Figure 12, the excitation is a node less than 4 inches from sensor D4. 
The proximity of the excitation to D4 is so small that accurate Time-of- Arrival 
estimates are difficult to determine. Locations D1 and D2 are equidistant from the 
excitation node with D3 significantly farther. The comparison of Time-of- Arrival 
versus geodesic distance is shown in Figure 13. In this case, the wave propagation 
speed is 2258 m/s. 

4. Summary 

Two global structural models of the on-orbit BEAM were generated and transient 
dynamic simulations of pseudo-MMOD impacts were completed. The model 
development supported understanding of MMOD impacts for Soft-Goods structures 
and specifically focused on restraint layer accelerations that would be measured by 
the on-orbit Distributed Impact Detection System (DIDS) to be deployed on BEAM. 
Considerable uncertainty exists as to the value for the elastic modulus of the 
restraint-layer. Fortunately, the parameters studies showed that the acceleration 
responses are nearly insensitive to the elastic modulus. However, the dispersive 
wave speeds for the acceleration responses are significantly dependent on the 
inflation pressures. The tension in the restraint-layer straps is directly related to the 
inflation pressure. 

Time histories at the restraint layer response locations were recorded for simulated 
impacts at randomly selected nodal locations on the Soft-Goods. Plots of Time-of- 
Arrival versus geodesic distance for multiple excitations indicated a wave speed of 
901 m/s. 

A number of modeling and measurement concerns have been identified for future 
effort: 

1. Incorporate a more realistic representation of the through-the-thickness 
Soft-Goods structure and impact physics utilizing results from Local BEAM 
Model studies, such as: 

a. Retain the current simple Soft-Goods Model and utilize a refined Local 
BEAM restraint-layer response. 

b. Map the complex transmission path of an MMOD impact through 
multiple layers and foam spacers. 
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2. Generate simulation data to multiple random impacts: 

a. Provide a library of response spanning a broad range of impact 
locations and severities. 

b. Provide sample signals to support DIDS data processing. 
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Figure 1. Graphic representation of BEAM attached to ISS. 
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Figure 2. NASTRAN Inflated Model. 


11 




(a] Soft-Goods Only 


[b] All Parts Except Soft-Goods 


Figure 3. Soft-Goods Model derived from Inflated Model. 
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Interior of Bulkhead Model 


Exterior of Bulkhead Model 


Figure 4. Schematic of Bulkhead Model. 
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Figure 5. Comparison of first 3 free vibration modes. 
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Figure 6. Effect of Inflation pressure on acceleration wave for Soft-Goods Model. 
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Figure 7. Effect of elastic modulus on acceleration response for Soft-Goods Model. 
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Figure 8. DIDS accelerometer locations. 
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Figure 9. DIDS acceleration locations along with nodal excitation locations shown in 

green. 
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Figure 10. Sample time histories for Soft-Goods impulse at the Soft-Goods DIDS 

sensor locations. 
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Figure 11. Correlation of impulse Time-of-Arrival with geodesic distance for Soft- 

Goods excitations and responses. 
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Figure 12. Sample time histories for bulkhead impulse at the DIDS bulkhead sensor 

locations. 
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Figure 13. Correlation of Time-of-Arrival with geodesic distance for bulkhead 

excitation and response. 
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